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A survey is given on the present status of the nuclcon parton distributions and re- 
lated precision calculations and precision measurements of the strong coupling constant 
as{M'^). We also discuss the impact of these quantities on precision observables at 
hadron colliders. 
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^\i . 1. Inside Nucleons 

> 

p\J ' The physics of the strong interactions always has been tightly connected to the 

C^ ■ study of nucleons at shorter and shorter distances. The measurement of the anoma- 

lous magnetic moments of the protorP and neutrorP in 1933 and 1939 made clear 
that nucleons are no elementary particles. During the 1950ies the Hofstadter exper- 

i— ^ . imenta^ revealed the charge distributions inside nucleons^ at scales Q^ ~ 0.5 • Mfj. 

f^ I Yet it was unknown how these distributions came about. In 1964 Murray Gell- 

ManrP proposed the quark model, to catalog the plethora of observed baryons and 
mesons. Independently G. Zweig suggested aces^ as the building blocks of hadrons. 
A direct connection to the lepton-nucleon scattering data was not made at that 
time. 
5^ \ Back in 1954 C.N. Yang and R Millg^ proposed novel bosonic field theories 

based on gauge invariance with respect to non-abelian groups. This development 
went unrelated to strong interactions for a long time. With the advent of the Stan- 
ford Linear Accelerator in 1968 the nucleon structure could be resolved at much 
shorter distances by the MIT-SLAC experimenta^"i^ beyond the resonant region 
W > 2GeV for values Q^ up to 30 GeV^. The remarkable finding by these experi- 
ments were that i) the structure function vW2{v,Q^) which has been expected to 
depend on both kinematic variables v and Q^ independently, turned out to take the 
same values for fixed values of x = Q"^ / {2M ^v) irrespectively of v and Q^ at high 
enough values. This phenomenon is called scaling, ii) The ratio of the longitudinal 
structure function Wl and W2 turned out to be very small. BjorkerP^ had predicted 



* Dedicated to M. Gell-Mann on the occasion of his 80th birthday. 



scaling at asymptotic scales Q^,v -^ oo in 1969. Learning about the SLAC-MIT 
results R. Feynman very quickly proposed the parton modet^, which is equivalent 
to Bjorken's description but based on the observed strict microscopic correlation 
between Q^ and v = q.pi 

W{., Q-) = E ^^ f dxM-^)5 (if - |i) , (1) 

where e^ and fi denote the parton's charge and distribution functions. Would the 
parton model be unique in describing the new data? This has been challenged by 
other popular formalisms like vector meson dominancei^. However, they failed to 
describe the behaviour observed for Wl, which corresponded to that of spin 1/2 
partons, according to the calculations by Callan and Gross^^. 

Yang-Mills theories^ became building blocks of the electro-weak Standard 
ModeP^, although there renormalizibility had not been proven yet, a conditio sine 
qua non for a physical theory. The proof was an urgent matter and in 1971 it was 
achieved both for massless and spontaneously broken Yang-Mills theories, along 
with designing practical loop computations in this sophisticated theory in an au- 
tomated way^SESI^ Quantum Chromodynamics (QCD) was proposed as the theory 
of the strong interactions in 1972 by Gell-Mann and Fritzsch and Leutwyler^ZISll 
as a renormalized Yang-Mills field theory based on S'C/(3) gauge interactiona^Sl. 
D. Gross, F. WilczelP^and D. Politzer^ studied the asymptotic behaviour of color 
octet gluon Yang-Mills theory, cf. alscP^I, and found asymptotic freedom. This is the 
essential ingredient, which makes it possible to perform perturbative calculations 
at large scales in a theory with strong interactions at low scales. 

At short distances the nucleon structure functions Fi{x, Q"^) obey the light-cone 
expansiorp31. At large scales Q^ the contributions of lowest twist dominate and the 
representation 

F,{x,Q^) = E^'(^'QV/^') ® /j(2^,m') (2) 

3 

holds. Here C| (x, Q"^ / fj?') denote the Wilson coefficients and fj [x, fj^) are the parton 
densities, /i^ is an arbitrary factorization scale and ® denotes the Mellin convolution. 
The scale behavior of the nucleon structure functions Fi{x^ Q^) obey renormal- 
ization group equations, an important aspect of renormalizable Quantum Field The- 
ories to which Murray Gell-Mann made very essential contributions very early^^. 
Transforming Eq. ^ to Mellin space one obtains the following Callan-SymanzilPSl 
equations : 



^|+^(^)|-'^^(^) 



F,[N,Q'')^Q (3) 



^It is interesting to note that different approaches to renormalization result into different mathe- 
matical structures as shown irl2SI,Thus the method by Gell-Mann and Lo»22Iis in general related 
to a cocycle, while that by Stiickelberg and PetermannSSl relates to a group. I thank A. Petermann 
for pointing out Ref.ufil to me. 






'^|l+^(^)|-^'^(^'^) 



MN,ti') = (4) 

C'=(iV,QVM') = 0, (5) 



with /3(g) the QCD /3- function and 7k(-/V, /i) the anomalous dimensions. Both func- 
tions imply the scaling violations of the structure functions^ With progressing time 
the measurement of the deep-inelastic structure functions improved considerably 
and after 40 years, e.g. the precision of the structure function F2'^^{x,Q^) reached 
1% over a very wide range, ci!^. Due to this both the precision measurement of 
the unpolarized parton distributions and the strong coupling constant as{M'^) is 
possible from these data. 

2. Higher Order QCD Corrections to Deep-Inelastic Scattering 

On the theory side, the progress in higher order computations, likewise, has been 
enormous during the same period. The initial 1-loop resultP^'^^'^ are now widely 
improved to 3-loop order and even somewhat beyond. This is necessary to comply 
with the current precision of data. The status of the theory of deep-inelastic scat- 
tering is illustrated in the flowchart below. Here the dates indicate the year in which 
the corresponding correction to the respective quantities has been calculated. 

Let me mention the most far reaching results. For the QCD /3-function the 4- 
loop corrections were first computed by Vermaseren et ahl^in 1997. The unpolarized 
anomalous dimensions and massless Wilson coefficients are known to 3-loops for a 
series of momcnts^^ and in complete forrrpZHSH since 2004/05. A first moment for 
the non-singlct+ anomalous dimensions has been computed at 4-loops^ in 2006 
and more moments are in preparation. The heavy flavor Wilson coefficients in the 
region Q^/m^ = p ^ 1, which is a good approximation in case of F2{x,Q'^) for 
yO > 10, were computed for a larger number of momenta^ and in complete form for 
F]^. Currently the computation of the Wilson coefficients at general values of A'' 
is underwaySH. 

At the level of the leading twist (r = 2) representation the light cone expansion 
and the QCD-improved parton model lead to the same results. In the 1980ies it was 
thought, that the higher order corrections need to be supplemented by different 
small- a: resummations, cf! ^^ * ^^ l, to obtain correct results, even in the HERA kine- 
matic region. These perturbative resummations are connected to the problem how, 
within their approach, perturbative and non-perturbative contributions are clearly 
separated - an important pre-requisite to apply perturbation theory at all. 



''We remind that Drell and collaborators at the end of the 1960ies^ were seeking desperately 
scaling in fermion-meson interactions with loop corrections but ended up with scaling violations 
in general. 
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The resummatiorP^Ml j^g^g successfully predicted the so-called 'leading' poles of the 
QCD anomalous dimensions related to the poles at TV = 1 in Mellin space, at least 
up to 0{af). The leading series is related to the scale-invariant limit of QCD. The 
corresponding resummed anomalous dimension, however, has branch cuts in the 
complex plane^^but no poles at all, see alsoSSl. These singularities are much milder. 
Phenomenological studies^ have shown, that subleading effects are as important 
as the leading ones, since they widely cancel the effect of the former. One estimates 
that about four complete series of these terms are needed to obtain convergence. 
Currently the only practical approach relies on the computation of the Wilson- 
coefficients and anomalous dimensions to high enough order, which includes all 
the small- and large-x effects automatically. In the latter case, the renormalization 
group even allows reliable resummation^^. 

Beyond the level of leading twist much less is known on deeply-inelastic structure 
functions. Most of the results obtained so far concern the 1-loop level, cf. e.g.'^^. 
Here, the corresponding partonic operator matrix elements depend on several di- 
mensionless invariants Xi, unlike in the case of lowest twist. They cannot be mea- 
sured individually in the deep-inelastic process, but require ab-initio determinations 
using reliable non-perturbative methods. For a series of moments, this may be pos- 
sible in the future, using lattice techniques. In the polarized case a series of results 
has been obtained for the twist-3 contributions, cf. e.g.^^, among them a series of 
integral relations between different structure function3^Sl8ll_ Finally, deep-inelastic 
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Fig. 1. Left panels: The parton densities xu^ and xd^ at the input scale Qg = 4.0 GeV^ (solid 
line) compared to results obtained from NNLO analyses by MRST04 (dashed-dotted line)2^ and 
A06 (dashed line), 22 The shaded areas represent the fully correlated \a statistical error bands. 
Right panels; Comparison of the same parton densities at different orders in QCD from LO to 
N^LO Ref.123. 

non-forward scattering has been extensively studied during the last two decadespSl. 
These methods do in principle allow to measure the (quark) angular momentum of 
nucleonsSS, which is still difficult experimentally. One may apply these techniques 
to describe inclusive deep-inelastic diffractive scattering, referring to the light cone 
expansioiP^l and proving that the anomalous dimensions are structurally the same 
as in the forward case. 
3. Precision Quark and Gluon Twist-2 Distributions 

In the following we describe recent extractions of the unpolarized twist-2 parton 
densities to 3-loop accuracy and higher. 



3.1. Flavor non-singlet analysis 

The flavor non-singlet parton distributions obey scalar evolution equations and do 
not depend on the gluon distribution, which is more difficult to access in deep- 
inelastic scattering and may cause some systematic uncertainty, in particular deter- 
mining the strong coupling constant as(M§). Moreover, in the small x region the 
QCD evolution leads to moderate changes of the the distributions, unlike in the fla- 
vor singlet and gluon-case. Due to this flavor non-singlet analyses are advantageous. 



One may apply the valence-approximation in the region x > xq, xq ~ 0.35, 0.4 and 
construct a non-singlet distribution from deuteron and proton data for a; < Xq. To 
describe the valence quark parton densities, the distribution x{d ~ u){x,Q'^) has 
to be known, which can be measured using Drell-Yan data. Furthermore, the non- 
singlet 0{al) heavy flavor corrections are applied, which amount to about 1%. The 
yet unknown 3-loop corrections will be even smaller. To perform a leading twist 
analysis, kinematic regions with higher twist contributions are cut out in a system- 
atic study, implying the cuts of Q^ > 4 GeY^,W'^ > 12.5 GeV^, dW^. The results 
for the parton distributions xUy{x,Q^) and xdy{x,Q'^) are illustrated in Figure 1 
and are compared to other determinations. Wile in the case of the xUtj distributions 
the overall agreement is good, there are still systematic differences in case of the 
down-valence distribution. We also illustrate the perturbative expansion from LO 
to N'^LO reaching convergence. 

Extrapolating the twist-2 QCD fit results into the region 12.5 GeV^ > W"^ > 
4 GeV^ the flavor non-singlet higher twist contributions can be determined empir- 
ically^H. The inclusion of soft resummation terms for the Wilson coefficient beyond 
the N'^LO corrections allows to extract the higher twist terms in the region x < 0.75 
in a stable way, while lower order analyses overestimate the higher twist contribu- 
tions. The relative higher twist contributions in the proton and deuteron case turn 
out to be about of the same size. 

3.2. Combined singlet and non-singlet analysis 

In combined singlet and non-singlet analyses of the deep-inelastic world data at 
NNLO, cf. Refs.'2S"22, one determines also the different sea-quark and gluon densi- 
ties. This has always to be done together with the measurement of the QCD scale 
Aqcd due to strong correlations. To unfold the sea-quark densities one refers to 
Drell-Yan- and di-muon data as well, through which the distributions x{d—u)(x, Q^) 
and xs{x,Q^) — xs{x,Q^) can be measured individually. Due to the large charm- 
quark contribution to the deep-inelastic structure functions the description of the 
heavy flavor contributions is required at the same level of accuracy as for the light 
partons. Currently it is available to 0{al'p^ and the 0{al) corrections are under- 
way^iESESSl^ To obtain an accurate interpolation between low and higher scales Q^, 
the so-called BMSN-interpolation is recommendecP^lEnil, Nowhere in the kinematic 
region of HERA heavy flavor logs become large to be resummecP^, i.e. in a very 
wide kinematic region even the charm quarks cannot be viewed massless. Yet, one 
may define heavy flavor parton densities in terms of technical quantitie j^^ * ^°^l to 
some extent 1^ to evaluate other observables. Here one has always to check to which 
extent this approximation holds. 

The present analyses have lead to very precise parton distributions. The NNLO 
parton distributions determined in Refs.'^SES do widely agree within the measured 



'^At 3-loop order graphs exist with both charm and bottom-quark lines in the operator matrix ele- 
ments. They do not fall under the paradigm of single parton distributions, despite being universal. 
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Fig. 2. The light parton densities xu, xd and xG at the scales /^^ = 4, 100, 10000 GeV^ . The bands 
denote the parton distributions with la uncertainty of ABKMOSJ^SJ, The dashed lines correspond 
to MSTW0&29.; from Ref.i^. 



region, with very slight differences in the x(u + u) and x{d + d) distributions. In 
Figure 2 we compare the results of the NNLO fits of Refs'^^'^for the light partons. 
At low scales fj,^ the sea-quark and gluon distributions^^ take lower values than 
those in Ref.'^Sl, and the NNLO gluon distribution even tends to negative values 
yielding the largest difference. 

Through the evolution the densities get closer. This, however, is partly due to 
the large value of as{M'^), Eq. ([14]), which leads to a relative acceleration of the 
evolution compared tcP^I. At scales larger than Q^ ~ 10* GeV^, accessible at the 
LHC, this will lead to a further growing gluon density o^ compared tcPS. The 
precision observables at the LHC will help to constrain the parton distributions 
further. 



4. The Strong Coupling Constant 

The QCD parameter Aqcd, or as{M'^), is determined in QCD fits together with 
the non-perturbative input densities for the different partons at a starting scale 
Qq. There are tight correlations between the value of as{M'^) and some parameters 
of the parton densities. An important example is the normalization of the gluon 
density, ciW^^^. In the non-singlet analysi^we obtained at NNLO, cf. alsJ^^Zl^ 
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226 ± 25 MeV 



(6) 



and at N'^LO, assigning to the yet unknown 4-loop anomalous dimension a ±100% 
error, 
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Usually the QCD parameter is expressed in terms of as(Af^)- In the following we 
compare the results of recent NNLO and N'^LO analyses for the deep- inelastic world 
data obtained by different groups : 

a.(M|) = 0.1134 +;!-';';i? NNLO m (8) 



-0.0021 
a.(M|) = 0.1141 +||;^|]20 ^3lo ^ (9) 

as{Ml) = 0.1135 ± 0.0014 NNLO, FFS 
a,{M^) = 0.1129 ± 0.0014 NNLO, BSMN 
as(Af|) = 0.1124 ±0.0020 NNLO, dyn. approach 
a,.(Af|) = 0.1158 ± 0.0035 NNLO, stand, approach 
a,.(Af|) = 0.1171 ± 0.0014 NNLO 
More recent unpolarized NNLO fits, including the combined HERA datcP^, yield 
as (Af|) =0.1147 ±0.0012 NNLO [HS] (15) 

asiM^) = 0.1145 ± 0.0042 NNLO, preliminary [109] (16) 

Note that the values (|8|9p are independent of the gluon distribution. The combined 
singlet non-singlet analysis, based on rather different data, and being sensitive to 
both the sea-quark and gluon densities, yield the very similar values (|10|lip . This 
analysis has been performed with a different code than used irF^l. 

The above values are located below the present weighted average of as(Af^) 
measurementsii^ of 

as(Af|) = 0.1184 ± 0.0007 , (17) 

cf. Figure 3. The error given in (|17p cannot include the yet unknown relative sys- 
tematics between the different classes of the same type of measurement. 

We would like to mention that recent determinations of as{M'^) using event 
shape moments for high energy e+e^ annihilation data from PETRA and LEP 
including power corrections the following values were obtained : 

as{Ml) = 0.1135 ± 0.0002 (exp) ± 0.005 (f^i) ± 0.0009 (pert) NNLO [111] (18) 
as{Ml) = 0.1153 ± 0.0017 (exp) ± 0.0023 (th) NNLO [Il2]. (19) 

Also these measurements of as{M'^) yield low values. They show that the results 
obtained analyzing deep-inelastic data do not form a special case. Also in deep- 
inelastic scattering off polarized targets as{M'^), at NLO, has been measured, how- 
ever with larger errors, sec^^^^^' ^^^-'^'^l The present error on as{M^) at NNLO of 
~ 0.0012 is at the margin of the present theory and systematics errors. Different 
known theoretical uncertainties, cf. also Eqs. (|8l-[TT]). are of the order of 0.0007, the 
quoted la error of the world average. The systematics of the different extractions 
of as(Aff ) has to be understood in even more detail in the future. The current val- 
ues of as(Af^) obtained from precision deep- inelastic scattering data disfavor the 
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Fig. 3. The strong coupling constant ciaiKi^) from diff erent PI S measurements, at NNLO and 
N^LO, Eqs. lIMTSll. and at NLO in the polarized case, cf.ESlESl. The yellow band describes the 
weighted average of a wide range of (is(Mz) measurementaH^J RefPSSl 

unification of forces, even in the supersymmetric extension of the Standard Model. 
However, the picture may change soon, with new findings at the LHC. 

5. Main Inclusive Cross Sections at Hadron Colliders 

The accuracy for the parton distribution functions reached can now be apphed to 
derive precision predictions for inclusive hadronic observables, such as the Drell-Yan 
cross section, the W^ , Z-boson, the ti- and Higgs-boson production cross sections 
at NNLO. Detailed analyses have been given in^'^^^^. For all these quantities at least 
these corrections are necessary. The Drell-Yan cross section and the W^ , Z-boson 
production cross sections are, furthermore, used as 'standard candle' processes to 
measure the collider luminosity. They have therefore to be known as precisely as 
possible. As an example we show in Figure 4 predictions on the inclusive Higgs 
boson production at hadron colliders. Within the present accuracy of the parton 
distribution functions the predictiongSSESI g^ju show some differences, which are 
likely related to the different gluon distributions at low scale and the values of 
as{M\). The predictions agree for s ~ Mz but differ for higher mass scales. 



6. Higher Loop Integrals and Mathematics 

The computation of higher order loop integrals is still a difficult task, even in the 
massless case or in the presence of a single mass scale. At three-loop (and higher) 
orders this both applies to the zero- and single-scale problems calculated at present. 
At the one hand, the results expected have a rather simple structure. On the other 
hand, a growing multitude of diagrams which contain more and more difficult struc- 
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Fig. 4. Inclusive Higgs boson production cross section at TEVATRON and the LHC. The bands 
denote the parton distribution uncertainty of ABKMOS^^; the lines correspond to MSTW0823 
from Ref.EI!. 



tures have to be computed. Obviously, an enormous part of intermediary results 
simply cancels. It is, however, difficult to let cancel these contributions at a rather 
early stage of the computation, or to even widely avoid that they occur from the 
very beginning. This is one of the central problems of all present calculations. Gauss' 
theorerrP31] allows to express Feynman diagrams to a set of master integrals, which 
finally have to be computed for zero scale quantities. Also in case of single scale 
quantities, e.g. given by a Mellin variable A'', one may obtain similar recursions. This 
method, simplifying the calculation, may lead to a large number of terms being of 
higher complexity than those finally appearing in the results. This property seems 
to be in common with different other approaches, which in the first place appear 
to simplify the calculation technically, like Mellin-Barnea^^ integrals or multino- 
mial expansion, since they lead to rather elementary Feynman parameter integrals. 
Another approach, cf. e.g. 2312211211122]^ consists in evaluating the individual Feyn- 
man parameter integrals without applying intermediary simplifying methods, i.e. 
mapping them directly to the analytic mathematical structure they represent. This 
method is more demanding but will finally lead us to a deeper understanding of 
the objects we deal with. For single scale single integrals integrals and one mass, 
the 2-loop integrals and simpler 3-loop integrals can be represented in terms of 
generalized hypergeometric functions and extensions thereof, like Appell functions 
(,£li23| j^^ important issue in integration is a clear definition of the target space and 
the knowledge of the relations of its elements. For the Feynman parameter inte- 
grals discussed above, shuffle- and Hopf-algebraa^^S play a central role, along with 
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Poincare iterated integrals over specific alphabets. Feynman integrals are almost 
always periods^^-. In this context, the simplest structures are multiple zeta val- 
ues, cfP^, nested harmonic suma^^, harmonic polylogarithm^i^il^ g^jjj^ generalized 
harmonic sum^iSSHlsT] Peynman parameter integrals will coin new classes of higher 
transcendental functions going to higher and higher order, which have to be studied 
to perform future precision calculations in an efficient way. Their evaluation is in- 
timately connected to modern summation technologies, like SIGMA^'^^, and efficient 
algorithms to establish and solve the associated recurrences of both large order and 
degre^^21i39| ^11 this will require high performance computer algebra written us- 
ing highly efficient languages like F0Rl4i^, nearly 50 years after SCHDDNSCHIP was 
introducecP^, and the investment of many CPU years, however, at an even more in- 
volved level than considered today. Structures of Feynman parameter integrals, e.g. 
on the level of multiple zeta values, form an interesting recent field in mathematics, 
^,£j i22 | i4i [i42j^ which is also related to irrationality proofs of the basis elements span- 
ning these quantities. Here we face a new era of a tight symbiosis between theoretical 
physics and modern mathematics, which is regarded to be very essential. 

7. Conclusions 

Quantum Chromodynamics was a great discovery. With it Murray Gell-Mann com- 
pleted the revolution of the strong interactions started in the early 1960ies with the 
introduction of the quarks. During the last 37 years computations grew to a pre- 
cision of 0(1%) for inclusive quantities, which are described at 3-loop, and partly 
at 4-loop, level, moving the frontiers of Quantum Field Theory to breathtaking 
new horizons. The running of the strong coupling constant is understood in great 
detail, despite in different classes of analyses still values of q;s(M|) are found which 
differ by experimental and theoretic systematic effects, being partly yet unknown. 
To determine as{M^) at the level of its present statistical accuracy of ~ 1% further 
studies and even higher order calculations are required for some of the processes. 
The QCD improved quark-gluon parton model works impressively well at short dis- 
tances - a clear triumph of Quantum Chromodynamics and proof that quarks and 
gluons, although being confined, are basic building blocks of matter. Without them 
the Standard Model would suffer from anomalies^^ and not form a Quantum Field 
Theory. 

During the last two decades the methods of lattice QCD steadily improved, in 
particular concerning the systematic errors involved. We therefore expect precision 
computations both on Aqcd and a series of moments of unpolarized and polarized 
parton densities in the near future. The results of these calculations ab initio can 
then be compared to the precision extractions discussed based on precision exper- 
imental data and higher order perturbative calculations. A final question concerns 
the unification of the three forces of the Standard Model, given what we know at low 
scales at present. Future discoveries, perhaps at the LHC, will lead to a clarification 
here. 
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